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Abstract 

This  is  a  report  describing  the  research  performed  during  the  past  two  years.  The 
project  has  been  funded  for  two  years,  and  additional  funding  for  a  third  year  is  requested. 
The  contents  of  this  report  includes  the  following:  (I)  a  description  of  the  typical  ultrasonic 
submersion  test  set-up,  (H)  a  brief  summary  concerning  the  research  achievements  for  the 
past  two  years.  A  renewal  of  one  year  was  requested  to  complete  work  already  brought  to 
a  high  level  during  the  past  two  years  and  is  currently  being  funded. 

I.  Typical  ultrasonic  submersion  test  set-up 

The  ultrasonic  submersion  test  is  a  useful  technique  for  the  inspection  of  material 
properties.  The  present  research  topic  concerns  the  use  of  the  ultrasonic  submersion  test 
as  a  means  to  inspect  multilayered  flber  reinforced  composite  plates.  Each  layer  of  the 
composite  is  a  transversely  isotropic  solid  characterized  by  flve  elastic  constants. 

The  signiflcance  of  an  ultrasonic  test  relies  on  its  ability  to  interpret  the  acoustical 
signals  picked  up  by  the  transducer.  Figure  ( 1)  illustrates  the  ultrasonic  submersion  test 
set-up  given  by  Wu  &  Ho  (1990).  A  transducer  sends  a  pressure  wave,  which  propagates 
through  the  coupling  fluid,  penetrates  thrraigh  a  sample  plate  and  teaches  the  receiver 
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transducer  located  on  the  opposite  side  of  the  coupling  fluid. 

The  goal  of  this  research  is  ultimately  to  determine  the  structural  properties  of  a 
sample  composite  plate,  such  as  its  layer  configuration  and  stiffness.  This  goal  can  be 
achieved  by  understanding  well  the  physical  principles  of  the  test  so  that  the  acoustical 
signals  picked  up  by  the  receiver  transducer  are  correctly  evaluated.  For  example,  the 
interpretation  of  the  arrival  time  of  acoustical  signals  is  related  to  the  path  of  a  wave 
package  and  its  speed,  both  of  which  depend  on  the  sample  stiffness.  The  interpretation  of 
the  pressure  level  in  the  receiver  transducer  is  related  to  the  transmission  coefficient, 
attenuation  and  geometrical  decay  of  the  propagating  waves.  The  transmission  coefficient 
depends  on  the  mechanical  properties  of  the  sample,  such  as  its  impedance  and  other 
factors.  The  geometrical  decay  is  related  to  the  propagating  distance  from  the  pitcher 
transducer  to  the  receiver  transducer.  A  successful  interpretation  of  an  ultrasonic  test 
means  that  such  stmctural  properties  of  the  sample  as  its  layer  configuration  or  stiffness 
may  be  precisely  extracted. 

Because  of  its  close  relation  with  the  stiffness  of  the  sample  and  its  easy 
interpretation  from  the  picked  up  signals  (wavefront),  the  determination  of  the  wavespeed 
is  one  of  the  major  concerns  in  an  ultrasonic  test .  Two  types  of  wavespeeds,  body-wave 
and  guided- wave,  are  the  main  concern  in  a  test. 

When  a  high  ficquencyyfave  package  penetrates  through  a  composite  plate,  the 
wavespeed  of  major  concern  is  the  body  wave  type.  The  first  distinctive  feature  of  waves 
in  anisotropic  media  is  the  difference  between  the  direction  of  power  flow  and  the  direction 
of  phase  flow.  In  figure  2  (Auld,  1973),  the  direction  of  phase  velocity  is  horizontal  while 
the  direction  of  power  flow  is  not.  In  anisotropic  media,  the  path  of  a  wave  train  is 
interpreted  to  mean  the  path  of  power  flow.  The  ray  path  in  the  composite  plate,  figure  1 , 
is  the  direction  of  power  flow,  not  the  direction  of  the  phase  velocity.  The  path  and  the 


2 


velocity  of  power  flow  can  be  easily  measured  from  the  location  of  the  receiver  transducer 
and  the  arrival  time  of  the  wavefront  of  a  wave  package.  Wu  and  Ho  (1990)  successfully 
interpreted  the  location  of  the  transducer  and  the  arrival  time  of  the  acoustical  signals  to 
obtain  the  flve  elastic  constants  of  a  single  layered  composite.  However,  they  did  not 
study  multilayered  plates  and  did  not  study  the  propagating  (decay)  feature  of  waves  in  the 
composite. 

When  a  low  frequency  wave  package  is  incident  from  certain  special  angle  and 
excites  a  guided  wave  mode,  the  wavespeed  of  concern  is  a  guided-wave  type  (Lamb  or 
leaky  Lamb  wave).  The  guided  wavespeed  is  obtained  from  solving  a  frequency  equation, 
which  is  an  equation  giving  the  relationship  between  wavenumbers  and  excited  frequency. 
For  multilayered  media,  the  frequency  equation  can  be  found  from  the  propagator  matrix  of 
the  system.  Many  recent  articles  published  in  Journal  of  the  Acoustical  Society  of  America 
fall  into  the  category  of  investigating  the  guided  wavespeed,  see  Chimenti  &  Nayfeh 
(1988,1989,1990),  Dayal  &  Kinia(1989),  Karim,  Mai  &  Bar-Cohen  (1990). 

The  main  difflculty  involved  with  interpreting  the  pressure  level  picked  up  by  a 
receiver  transducer  is  taking  into  account  the  effect  of  attenuation  due  to  imperfection  and 
geometrical  decay.  As  the  slowness  surface  for  waves  in  anisotropic  media  may  not  be  in 
a  convex  shape,  a  focussing  effect  may  be  induced  in  the  wavefront  surface.  Thus,  the 
second  distinct  feature  of  waves  in  anisotropic  media  is  that  the  geometrical  decay  of  waves 
may  not  be  1/R  in  a  three-dimensional  space  and  1/R^  in  a  two  dimensional  space. 
Without  excluding  the  effect  of  geometrical  decay  from  the  measured  signals,  a  wrong 
conclusion  may  be  reached  concerning  the  attenuation  caused  by  the  imperfection. 

Because  a  composite  plate  mainly  serves  the  structure  by  its  use  as  a  plate  or  a 
membrane,  it  is  also  of  interest  to  connect  the  guided  (Lamb)  wave  modes  to  the  strucniral 
vibration  modes.  The  wavespeed  in  structural  vibration  modes  can  be  analytically  obtained 
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and  will  provide  the  asymptotic  solutions  for  the  dispersion  equation. 

In  summary,  the  essence  of  the  ultrasonic  testing  is  the  ability  to  interpret  the 
received  acoustical  signals.  However,  this  ability  relies  on  the  understanding  of  the 
physical  principles  about  waves  passing  through  the  testing  sample.  Thus,  wavefront 
analysis,  reflection  and  transmission  coeflicients,  guided  waves  and  geometrioil  decay  are 
major  topics  in  this  research.  The  achievements  made  in  understanding  those  topics  during 
the  past  two  years  of  research  will  be  discussed  next. 


II.  Summary  of  Current  Achievements 
(A)  Guided  waves. 

The  first  part  of  the  result  relates  to  the  understanding  of  guided  waves.  The  far- 
fleld  displacement  expressed  in  terms  of  double  Fourier  integrals  is: 


Uk 
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fk(n  i.n2>fa)) 
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exp[i  (TllXi+Tl2X2-(Ot)]  dnidTi2 


(1) 


The  poles  of  the  Fourier  integrals  are  related  to  guided  waves  and  g(r)  j  =  0 
is  called  the  frequency  equation.  The  functions  fj^  are  due  to  excitations  or  loadings  in  a 


small  region.  The  following  represents  the  conclusions  of  the  past  two  years; 

(1)  The  dispersion  function  g  can  be  constructed  by  propagator  matrix  method. 

(2)  For  a  single  layer  plate,  there  are  several  asymptotic  solutions  for  the  frequency 
equation  g  ■■  0.  At  low  frequencies,  three  asymptotic  solutions  are  constructed  by 
membrane  and  KirchhofT-plate  models.  At  high  frequencies,  three  asymptotic  solutions  are 
obtained  from  the  plane-strain  model  and  the  SV  wave  model.  A  special  solution  is  related 
to  the  common  point  of  the  slowness  surfaces  of  quasi-SV-wave  branch  and  SH-wave 
branch. 


(3)  The  wave  forms  at  the  far  field  can  be  obtained  by  applying  method  of  residues 
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and  method  of  stationary  phase,  and  the  corresponding  result  is 


fk(ni.n2)cxp{  i  [niXi^-n2X2-(0t-«-sgn(K)7i/4-H7r/2]} 

(27r|ic|r)‘/2|Vg(ni,n2)l 


(2a) 


where  k  is  the  curvature  of  the  slowness  curve  g  =  0  for  a  specific  frequency.  Explicitly, 


anf  5ti2  omani  an,  an2  dr\l  ^n, 


(2b) 


The  expression  for  the  far- field  displacement  field  breaks  down  when  k  »  0,  which  occurs 


when  the  slowness  curve  does  not  have  a  convex  shape,  and  the  corresponded  guided 
wave  decays  slower  than 


t-x 


(4)  The  wavefront  curve  of  Lamb  waves  at  time  t  due  to  a  point  source  at 
I  -  X2  “  0  is  the  curve  Xj^  “  Vj^t  where  the  velocity  vector  of  the  power  flow  (Lamb 


wave)  is 

-  ag/ank 
a»0w 

(5)  One  limiting  case  for  which  an  analytical  formula  can  be  obtained  is  that  of  a 
layer  sufficiently  thin  that  the  layer  is  essentially  in  a  plane-stress  state.  For  a  thin  plate 
that  is  subjected  to  a  horizontal  (in-plane)  point  force,  the  equation  of  motion  is 

OijklGkni.lj+Sim^P  Gim 

where  Gj^  is  the  displacement  along  the  k  direction  due  to  a  point  load  applied  along  the  m 
direction  at  the  origin,  and 

stress  state.  The  steady-state  displacement  field  written  in  the  form  of  double  Fourier 
integrals  is 

f  (n,x,^n2xru>01  dn.dn, 

j.oo  * 


^ijkl"%1d*'iJ33'id3/‘3333  “ 


where  the  explicit  form  of  the  dispersion  function  is 
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2  22  2  22  2^2 
g  =  (qiiiini+qi2i2n2*P<^  )(Qi2i2ni+q2222n2'P<^  )-(qi2i2+qii22)  nin2 

From  g  =  0,  the  explicit  relationship  between  wavenumbers  and  frequency  is 
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where  c  =cos  0,  s  «sin  0  and  q  j=qcos0,  q2~nsin0.  Also,  dg/dx].  and  6  g/dq.dq.  are 
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oqi 
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(B)  Body  Waves: 

The  far-field  displacement  expressed  in  terms  of  triple  Fourier  integrals  is: 

00 

Ffc(ni>n2.n3  )  (njXi+q2X2+Ti3“x3-wt)]  driidnzdTi 

CO  S“(qi,q2,n3) 


a 
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(3) 


For  this  case  the  major  concern  is  to  determine  how  the  incident  wave  package  propagates 
from  the  first  fluid-solid  interface  to  the  second  fluid-solid  interface.  The  propagation  of 
waves  from  the  second  fluid-solid  interface  to  the  receiver  transducer  is  not  of  immediate 
concern  as  the  waves  are  propagating  in  the  coupling  fluid.  The  function  F^  is  related  to 

excitations  or  loadings  in  a  small  region.  The  function  has  the  property  that  =  0  is 
the  equation  for  the  slowness  surfaces.  The  conclusions  of  the  study  are  given  next. 
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(1)  The  function  F®,  related  to  the  incident  waves,  can  be  obtained.  Carrying  out 
one  Fourier  integral  by  the  method  of  residues,  the  displacement  field  of  eq  (3)  becomes: 

mCO  mOO 

=  Lm  I  I  ?I^^?_expIi(niXi+n2X2+n3“x3-cot)]  dmdnidns* 

R—oo  J.oo  J.oo  ^ 

an? 

For  an  incident  wave  that  originates  from  the  fluid  side  and  penetrates  to  an  anisotropic  half 
space,  the  incident  wave  can  be  viewed  as  a  synthesis  of  plane  waves  with  amplitude 

A(npn2)- 


27riFg 

as“/an? 

-a. . 


=AT“d2  or  F?  = 


AT«dg  aS° 

27ri  an? 


where  T“  is  the  transmission  coefficient  of  a  plane  wave  and  d®  is  its  polarization  vector. 


(2)  The  wave  forms  of  eq  (3)  at  the  far  field  are 


iK«|‘^R  |VS«| 
where  the  Gaussian  curvature  is 


47r^H^ 

- - exp[i(Y+niXi+n2X2+n3"x3-<^0] 


[a^nfVa^f] 

\  ^ni  A  ^n2 , 

wni5n2/ 

fenff 

2 

and  Y  *  7t/2,  0  or  ir,  depending  on  the  signs  of  the  two  principal  curvatures.  The  point  of 
stationary  phase  (n  |  ,n2>n^)  satisfies 


S®(ni.n2.n3®)”0  and 
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X2 

X3 
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This  result  ako  means  that  when  the  phase  velocity  is  along  the  direction  (r|  j  ,ri2tn  ^ 
corresponding  direction  of  energy  velocity,  (Xj  ,X2,X2),  is  along  the  normal  direction  to  the 
slowness  surface,  i.e.,  (dS®/0rip  dS^/dr\2,  )• 

The  expression  giving  the  far  Held  displacement  breaks  down  when  curvature 
(k  -»  0).  This  case  will  occur  when  the  shape  of  a  slowness  surface  is  not  convex,  and  the 
decay  of  waves  is  slower  than  1/R.  Two  special  cases  are  related  to  cuspidal  edges,  which 
decay  as  or  conical  points,  which  decay  as  1/R*'^.  Our  study  confirms  that  the 

quasi-SV  slowness  surface  of  a  typical  flber  reinforced  composite  does  not  have  a  convex 
shape. 

(3)  Explicit  formulae  have  been  given  to  evaluate  the  required  flight  time  and  the 
horizontal  shift  of  the  receiver  transducer  to  pick  up  the  rays  for  a  multilayered  composite 
plate. 

(C)  The  reflection  and  transmission  coeflicients 

The  reflection  and  transmission  coefficients  for  a  plane  wave  incident  from  an 
arbitrary  angle  of  a  pitcher  transducer  to  the  receiver  transducer  were  also  studied.  The 
reflection  and  transmission  coefficients  were  obtained  by  the  method  of  propagator  matrix. 
The  effect  of  interfacial  debonding  is  also  considered.  Detailed  information  is  given  in  the 
paper  sent  out  to  publication  and  included  with  this  proposal. 

Publications  during  project  period 

1 .  W.  Lin  and  L.  M.  Keer,  “A  Study  of  Lamb  Waves  in  Anisotropic  plates”.  Journal  of 

the  Acoustical  Society  of  American,  in  press. 

2.  W.  Lin  and  L.  M.  Keer,  "Propagation  of  Waves  in  Fibre-Reinforced  Composite 

Plates”,  under  review. 


8 


References 

Auld,  B.A.  (1973)  Acoustic  Fields  and  Waves  in  Solids.  Wiley  Interscience,  New  York. 

V.T.,  Buchwald  (1959)  “Elastic  waves  in  anisotropic  media,”  Proceedings  of  Royal 
Society,  A253,  pp.563-580 

Chimenti,  D.E.  &  Nayfeh,  A.H.  (1989)  “Ultrasonic  leaky  waves  in  a  solid  plate  separating 
a  fluid  and  vacuum,”  J.  Acoust.  Soc.  Am.,  85-2,  pp.  555-560. 

Chimenti,  D.E.  &  Nayfeh,  A.H.  (1990)  “Ultrasonic  reflection  and  guided  wave 
propagation  in  biaxially  laminated  composite  plate,”  J.  Acoust.  Soc.  Am.,  87-4,  pp.  1409- 
1415. 

Cook,  T.  S.  and  Erdogan,  F.  (1972)  “Stresses  in  Bonded  Materials  wi.'i  a  Crack 
Perpendicular  to  the  Interface”,  Int.  J.  Engng.  Sci.,  10,  pp.677-697. 

Daniel,  I.  M.,  Anastassopoulos,  G.  and  Lee,  J-W.  (1989)  “Experimental  Micromechtmics 
of  Brittle-Matrix  Composites”,  in  Micromechanics:  Experimental  Techniques,  ed.  Sharpe 
Jr.,  W.  N.,  AMD  102,  ASME,  pp.  133-146. 

Dayal,  V.  and  Kinra,  V.  K.  (1989)  “Leaky  Lamb  waves  in  an  anisotropic  plate.  I:  an  exact 
solution  and  experiments,”  J.  Acoust.  Soc.  Am.,  85-6,  pp.2268-2276. 

Gao,  Y-C.,  Mai,  Y-W.,  and  Cotterell,  B.  (1988)  “Fracture  of  Fiber-Reinforced  Materials”, 
ZAMP,  39,pp.  550-572. 

Hsueh,  C-H.  ( 1988)  “Analytical  Evaluation  of  Interfacial  Shear  Strength  for  Fiber- 
Reinforced  Ceramic  composites”,  J.  Am.  Ceram.  Soc.,  7 1 ,  pp.  490-493. 

Hutchinson,  J.  W.  and  Jensen,  H.  M.  1990)  “Models  of  Fiber  debonding  and  Pullout  in 
Brittle  Composites  with  Friction”,  Mech.  Materials,  9,  pp.  139-163. 

Karim,  M.R.,  Mai,  A.K.  and  Bar-Cohen,  Y.  (1990)  “Inversion  of  leaky  wave  data  by 
simplex  algorithm,”  J.  Acoust.  Soc.  Am.,  88-  i ,  pp.482-49 1 . 

Marshall,  D.  B.  and  Evans,  A.  G.  (1985)  “Failure  Mechanisms  in  Ceramic  -  Fibei/ 
Ceramic-Matrix  Composites”,  J.  Am.  Ceramic  Soc.,  68,  pp.  225-231. 

McCartney,  L.  N.  (1989)  “New  Theoretical  Model  of  Stress  Transfer  between  Fibre  and 
Matrix  in  a  Uniaxially  Fibre-Reinforced  Composite”,  Proc.  R.  Soc.  Lond.,  A425,  pp. 
214-244. 

Nayfeh,  A.H.  &  Chimenti,  D.E.,  (1988)  “Propagation  of  guided  waves  in  fluid-coupled 
plates  of  fibre-reinforced  composite,”!.  Acoust.  Soc.  Am.,  83-5,  pp.  1736-1743. 

Sigl,  L.  S.  and  Evans,  A.  G.  (1989)  “Effects  of  Residual  Stress  and  Frictional  Sliding  on 
Cracking  and  Pull-Out  in  Brittle  Matrix  Composites”,  Mech.  Materials,  8,  pp.  1-12. 

Wijeycwickrcma,  A.  C.,  Kecr,  L.  M.,  Hiiashima,  K.  and  Mura,  T.  (1991a)  “The  Annular 
Crack  Surrounding  an  Elastic  Fiber  in  a  Tension  Field”,  Int.  J.  Solids  Structures,  27 ,  pp. 


9 


315-328. 


Wijeyewickrema,  A.  C.  and  Keer,  L.  M.  (1991b)  “Matrix  Fracture  in  Brittle  Matrix  Fiber- 
Reinforced  Composites”,  Int.  J.  Solids  Structures,  28,  pp.  43-65. 

Wijeyewickrema,  A.  C.  and  Keer,  L.  M.  ( 199  Ic)  “Matrix  Crack  Interaction  in  a  Brittle 
Matrix  Fiber-Reinforced  Composite”,  Int.  J.  Solids  Structures,  in  press. 

Wu,  T.-T.  and  Ho,  Z.-H.  (1990)  “Anisotropic  wave  propagation  and  its  application  to 
NDE  of  composite  materials,”  Experimental  Mechanics,  pp.313-318. 


Personnel 


During  the  project  period,  a  post-doctoral  fellow,  W.  C.  Lin,  assisted  L.  M.  Keer 
in  the  performance  of  the  research. 


10 


